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Abstract 

There has recently been considerable interest in the use of Nuclear Mag- 
netic Resonance (NMR) as a technology for the implementation of small 
quantum computers. These computers operate by the laws of quantum 
mechanics, rather than classical mechanics, and can be used to implement 
new quantum algorithms. Here we demonstrate how NMR can in principle 
be used to implement all the elements required to build quantum comput- 
ers, and briefly discuss the potential applications of insights from quantum 
logic to the development of novel pulse sequences with applications in more 
conventional NMR experiments. 

Key words: NMR, quantum computer, qubit, logic gate, controlled 
NOT. 
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1 Introduction 



It is well known that it is difficult to simulate the behaviour of a quantum me- 
chanical system with a classical computer. The difficulty arises because quantum 
systems are not confined to their eigenstates but can in general exist in any 
superposition of them; thus the vector space needed to describe the system is 
extremely large. For example a spin system comprising spin- 3 nuclei occupies 
a Hilbert space with 2^ dimensions. For this reason it is impractical to simulate 
the behaviour of spin systems with more than about a dozen nuclei. 

In 1982 Feynman ([J) reversed this observation, suggesting that quantum me- 
chanical systems have a potentially very large information processing ability. 
Thus it should be possible to build quantum mechanical computers which utilise 
this ability to achieve a computing power well beyond that of corresponding 
classical systems. The theory of such quantum computers is now fairly well un- 
derstood, but it has proved extremely difficult to actually build one. Recently, 
however, attempts to build computers based on the NMR properties of small 
molecules have exhibited considerable success (H, Y\, \T^, [71) • 

In this paper we show how NMR can be used to implement all the "compo- 
nents" required to construct quantum computers, and draw comparisons between 
quantum logic gates and more conventional NMR pulse sequences. The potential 
role of quantum computing as a source of new insights into NMR is also briefiy 
addressed. 

1.1 Bits and qubits 

The basic unit of information in a classical computer is the bit, which can take 
one of two values, and 1. Bits are then connected together by logic gates to 
form logic circuits, which can implement more complex logic operations, such as 
addition. Developments in classical computers have been driven by developments 
in the design and construction of logic gates, which have steadily become smaller, 
faster, and cheaper. However this process is beginning to reach a fundamental 
limit, as logic gates are reduced in size to atomic dimensions. Progress beyond 
this limit will require a different approach. 

One obvious possibility is to implement bits and logic gates using atomic 
components. A bit can be implemented using any two state device, such as the 
two quantum states of a two level system. For example the two Zeeman levels, 
\a) and |/?), of a nucleus in a magnetic field can be naturally described as a 
bit. Similarly a spin-system containing N distinct nuclei can be modelled as 
a set of bits. Traditionally the lower and upper energy levels are referred to 
as |0) and |1) respectively. 

The time-evolution of a spin-system system under some Hamiltonian is de- 
scribed by a series of unitary transformations, and so is of necessity reversible. 
Hence any quantum mechanical computer can only implement reversible opera- 
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tions, and must be built from reversible logic gates. This is not an important 
restriction as it has been shown that reversible logic gates can be used to simu- 
late traditional gates, and thus reversible computers are just as powerful as their 
irreversible counterparts 
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There is, however, much more to quantum computers than the implementation 
of classical algorithms using reversible logic: quantum computers are also capable 
of implementing new types of quantum mechanical algorithms ([I3,||^|7^|i^, with 
potentially enormous powers. This occurs because a two-level quantum system 
is not confined to its two eigenstates, but can exist in superpositions of these two 
states; that is the system is not confined to |0) and but can exist in states 
such as 



Co|0) + Ci|l) [1] 

where Cq and Ci are complex numbers and CqCq + c^Ci = 1. A nucleus in such a 
state is not really in state or state 1, but is in both states simultaneously. For 
this reason a two level quantum system is more than a simple bit, and is better 
described as a quantum mechanical bit, or qubit. A spin-system with N nuclei 
contains qubits, and can be in a superposition of up to 2^ states. This ability 
to be in a large number of states simultaneously gives quantum computers an 
intrinsic parallelism, which is exploited in quantum algorithms. 



1.2 Qubits and NMR spin states 

Traditional designs for quantum computers comprise A^ two-level systems which 
are coupled to one another and have some specific interaction with the outside 
world (so that they can be monitored and controlled) but are otherwise isolated. 
NMR systems are, by contrast, rather different. In particular a typical NMR 
sample comprises not just one spin-system, but a very large number of copies, 
one from each molecule in the sample. Thus while quantum computers are usually 
described using Dirac's bra and ket notation, NMR systems are described using 



a density matrix, usually written in the product operator basis ([7^). While it 
is possible to draw close analogies between the states of traditional quantum 
computers and those of NMR systems, it is necessary to proceed with some 
caution. 



1.2.1 One qubit states 

A single qubit can be in either of its two eigenstates, |0) and or in some linear 
superposition of them. Such a state is most conveniently written as a column 
vector in Hilbert space: for example the state described in equation |1] is written 
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as 

The corresponding density matrix 

can be decomposed as a sum of the four Pauh basis states, \E, I^, ly, and 

Consider first the eigenstates, |0) and These correspond to the density 
matrices 

|0)(0| = (J o)=|^ + /. [4] 

and |1)(1| = — respectively. Multiples of the unit matrix can be added to 
density matrices at will without effecting any NMR observable in any way, and 
so as far as any NMR experiment is concerned the density matrix 1^ is equivalent 
to |0), while —Iz is equivalent to |1). This simple approach is not, however, 
applicable to larger spin systems. 

Next consider superpositions, such as (|0) + |l))/\/2, with its corresponding 
density matrix 

(J + 

Once again multiples of the unit matrix can be ignored, and so is equivalent 
to (|0) + \1))/V2. Similarly |0) + is equivalent to ly, while |0) - |1) is 
equivalent to —Ix- Just as qubit eigenstates are closely related to magnetizations, 
superpositions are closely related to NMR coherences. 

1.2.2 Two QUBIT STATES 

While the relationship between qubit states and NMR states is simple for one 
qubit (one spin systems), this relationship is much more complicated in systems 
with two or more qubits. Indeed the problem of creating NMR states correspond- 
ing to multi-qubit eigenstates prevented progress in the implementation of NMR 
quantum computers for many years. 

Typically quantum algorithms start with all qubits in state |0), which for a 
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two-qubit computer is the state |00). The corresponding density matrix 



|00)(00| 



10 









I) 



[6] 



is quite different from the thermal equihbrium density matrix 
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-1 



) 



[7] 



Cory et al. have shown how this problem can be overcome (H,0, 0)- The ideal 



and this sum (ignoring multiples of the unit matrix as usual) can be assembled 
using conventional NMR techniques. An alternative approach, due to Gershenfeld 
and Chuang ([^), works by selecting four states from within a set of spin states 
arising from a larger spin-system. With a careful choice of states it is possible 
to find four levels whose relative populations correspond to those of |00)(00|, 
and these levels can be used as a pseudo two-spin system. While this approach 
is elegant, it is difficult to apply in practice and has not been widely used. A 
third approach, called temporal averaging ([^), is conceptually related to Cory's 
approach, but uses phase cycling instead of field gradients to select the desired 
state. 

Superpositions can be treated in much the same way, but they are not directly 
related to coherences in any very simple way. For example consider the state 
(|00) + |01))/\/2, in which the first spin is in state |0), while the second spin is 
in a superposition of states. The corresponding density matrix can be directly 
decomposed: 



but a more subtle approach is to note that (|00) + |01))/v^ can be written as a 
product of one-qubit states 



The corresponding density matrix can also be decomposed as a direct product of 





|00) + |01) _ |0)(|0) + |1)) 
V2 V2 



[10] 
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equations || and |^: 



'I 0' 
.0 0, 



1 1 

2 2 

1 1 

2 2 



{\E + /,) X {\E + S:) = \ [\E + h + S, + 2hS,) . [11] 



Note that, unlike the one qubit simple superposition does not correspond 

directly to an NMR coherence, but rather to a complex mixture of coherences 
and populations. Fortunately it is rarely necessary to directly consider issues of 
this kind, as such states can be easily obtained from states like Eq. |^. 

Finally we consider superpositions of the form (|00) + |ll))/v^, which cannot 
be broken down into a product of one qubit states (such states are normally 
referred to as entangled states). As such states cannot be factored it is necessary 
to decompose the corresponding density matrices directly. In this case 
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I ~ ^ (^"^ ~'~ ^-^-2*^2 ~^ '^■^xSx ^ '^■^ySy) 1 [12] 



i i 



which is a mixture of longitudinal two-spin order and DQ^ double quantum 
coherence. 



1.3 Global phase shifts 

One important consequence of the density matrix description of NMR quantum 
computers is that it is completely insensitive to global phase shifts. In general 
the wavefunction of any isolated system can be multiplied by an arbitrary phase 
shift without any observable consequences, that is the states and \ip') = e^'^lip) 
are indistinguishable. Indeed the absolute value of is completely meaningless, 
although it is possible to determine relative values of for otherwise indistin- 
guishable states by interference experiments. 

In the density matrix description of a state such global phases are not pre- 
served, as 

|^')(^'|=e^<^|^)(^|e-^^ = |V^)(^|. [13] 

Thus global phases have no discernable effect in any NMR experiment, and can 
be completely ignored. This is fortunate, as most NMR pulse sequences create 
global phase shifts, as discussed below, but as such phase shifts are truly global 
they can be neglected. 



2 One qubit gates 

One qubit gates act to modify the spin state of a single nucleus, and thus cor- 
respond to rotations in single-spin subspaces. Any rotation of this kind can be 
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achieved using RF fields, and so these gates are relatively straightforward. The 
gates can be implemented using selective pulses, in which case the gate is applied 
to a single nucleus, or using hard pulses, in which case the gate is simultaneously 
applied to a large number of separate nuclei. In this latter case the gate is more 
properly considered as a product of one qubit gates, one for each nucleus affected. 



2.1 The not gate 

The simplest one qubit gate is the not gate ([I^), which is well known from 
classical computing (thus not is a classical one bit gate, as well as a one qubit 
gate). This gate, which we shall call N, implements the rotation 

10) ^ ID 

|l>^|0) 

This can be described more compactly using a transformation matrix 

1 



where the symbol on the left signifies a not gate in a quantum circuit 
Clearly this gate can be implemented using a 180° Ja: pulse, as 

e-'"' = (_°, [16] 

and this matrix has the correct form up to a global phase change. This global 
phase term is irrelevant, as the overall phase is not an NMR observable quantity, 
and thus a 180°/a; pulse provides a good implementation of a not gate. 

As the NOT operation is simply an inversion, it is hardly surprising that it is 
implemented by a 180° inversion pulse. It might seem that any inversion pulse, 
such a 180°/j^, would be suitable, but this is not the case. For example 

performs the transformation 

— 1181 

ID ^ -|o> l^'l 

This is not a simple inversion, as it also negates the sign of the |1) state. This 
is not very important when the gate is applied to a system in an eigenstate, but 
is important when the gate is applied to a superposition. Consider, for example. 
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the superposition (|0) + |l))/-\/2, for which 

|0) + |1)^|0) + |1), [19] 

while a 1S0° ly pulse would give — 10) + ]1), a quite different state. This should 
not be surprising, as superpositions are closely related to NMR coherences, and 
Ix and ly pulses can have quite different effects, depending on the relative phases 
of the pulse and the coherent state. 

2.2 The square- root of not 

The square-root of not gate, V, is a purely quantum mechanical gate, in that it 
has no classical equivalent. As the name implies, V has the property 

= VV = N [20] 

and so an obvious implementation is a 90°/a; pulse, 

^""'■' = ^,{-i ?)• 

Once again this is equal to the "ideal" form, 

up to a global phase. 

The effect of V is to take an eigenstate to a superposition of eigenstates. For 
example 

|0)^(|0)-i|l))/V2. [23] 

This emphasises the quantum mechanical nature of V, as such superpositions do 
not have classical equivalents. 

Pulses with other flip angles can be treated in much the same way: for example 

a 60°Ix pulse is equivalent to a cube-root of NOT gate. This is not particularly 
interesting for one qubit gates, but becomes more interesting when comparing 
two qubit logic gates with spin state selective excitation sequences. 

2.3 The Hadamard gate 

The square root of NOT is not unique in converting eigenstates to superposi- 
tions: any 90° pulse will have a similar effect, as will a number of other pulse 
sequences. One particularly interesting sequence would be one corresponding to 
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the Hadamard gate, H, which performs the rotation 

|0)^(|0) + |1)) 



[24] 

|l)^(|0)-|l))/v^ 

This has two useful properties. First, it takes |0) to a completely uniform su- 
perposition, that is to a state where the coefficients in front of |0) and |1) are 
identical. Secondly, H is self inverse, so applying H twice is equivalent to doing 
nothing. 

H can be implemented in NMR using an off-resonance pulse, as 

)/v^ = _L 7), [25] 

which has the desired form, ignoring a global phase shift as usual. Alternatively 
this can be implemented using a three pulse sandwich ([^) , such as 

45 ly - 180 4 - 45 1.y. [26] 

When implementing quantum algorithms on NMR quantum computers, it is 
often easier to replace the Hadamard, H, by the pseudo-Hadamard operator, h, 
which has the form 



[27] 



This operation takes |0) to a uniform superposition of |0) and just like H, 
but it is not self inverse. This can be implemented using a 90° pulse. In many 
algorithms, a pair of H gates can be replaced by one h gate and one gate; this 
last gate is easily implemented as a 90° ^ pulse. Considered as an NMR operator, 
the H sequence performs the rotation 

^ 4 ^ /. [28] 
while h performs the more conventional rotation 

h^L^ -h [29] 



3 Two QUBIT GATES 

All the gates described above are one qubit gates: in NMR terms they perform 
rotations within the subspace corresponding to a single spin. If implemented 
using hard pulses they can in general effect several spins, but the overall effect 
is that each spin rotates within its own subspace, and the operation can be 
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decomposed as the product of two or more one spin operations. True two qubit 
gates correspond to rotations within a subspacc corresponding to two spins, and 
cannot be decomposed into a set of one qubit gates. These gates he at the heart 
of quantum computation, as they provide conditional dynamics: that is, the state 
of one spin can become dependent on the state of another spin. 

3.1 The controlled-not gate 

The fundamental two qubit gate is the controUed-NOT gate, which applies a not 
gate to one qubit (the "target") if another qubit (the "control") is in state |0). 
It can be described as follows: 



N 



'1 0^ 

10 

1 

,0 1 0> 



[30] 



A variety of methods for implementing this gate have been described, but it is 
more useful to consider a general approach to gates of this kind. It is well known 
in quantum computation that controlled gates, such as the controlled NOT, are 
related to controlled phase shifts by the Hadamard transform. For example a 
controlled- NOT gate can be replaced by the three gate network 



- N - 



- H - TT - H - 



[31] 



where H are one qubit Hadamard gates, and tt is given by 



TT 



a 

10 

10 

,0 -1, 



[32] 



which performs the transformation |11) while leaving other basis states 

unchanged. 

When implementing controlled gates in NMR pulse sequences, it is simpler 
to replace the Hadamard gates by pseudo Hadamard gates, as described above. 
Hence a general controlled gate can be implemented in NMR using the network 



-h-i- - h - 



[33] 
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where is a general controlled phase shift 




10 
1 
1 




[34] 







which performs |11) 6**^111), while leaving other basis states unchanged. 

This general form for controlled gates in NMR should not be a surprise, as 
there is a close link with composite 2;-pulses. It is well known that 2;-pulses can 
be replaced by three pulse sandwiches ([^, for example a 6^ pulse can be 
replaced by 90°2,6'y90°. Similarly, by cyclic permutation of 6x pulse can 

be replaced by 90°y6'290°, which is equivalent to \\~^6z^. Since a not gate (a 
180° pulse) can be implemented as an inverse pseudo Hadamard, followed by a 
180° phase shift, followed by a pseudo Hadamard, it is hardly surprising that a 
controlled-NOT gate can be implemented in much the same way, but using a 
controlled phase shift. 

3.2 Controlled phase shifts 

Controlled phase shifts, such as (f) (Eq. ^4|), are relatively simple to implement, as 
they can always be decomposed as a product of diagonal operators. For example 



The last three terms are straightforward, but the first term is difficult to obtain as 
it requires a Hamiltonian proportional to \E. This is not, however, important, as 
this term simply imposes a global phase shift, and as such can be ignored. For the 
remaining three terms, 2IzSz is proportional to the scalar coupling Hamiltonian, 
while Iz and Sz can be implemented as periods of free precession or by using 
composite z-pulses. 

The matrix tt, which lies at the heart of the controlled-NOT gate, can be 
implemented as 



which can itself be achieved in a variety of ways. The three terms above commute, 
and so the three Hamiltonians can be applied in any order, while the Iz and Sz 
terms can be implemented by free precession or by any of a wide variety of 
different composite pulses. 

Just like simple one qubit gates, two qubit controlled gates can also introduce 
global phase shifts, but as long as these are global and universal, that is they 
are applied to the the whole wavefunction and not just to the spins participating 
in the gate and they are applied irrespective of the state of the control bit, such 
phase shifts can be ignored. This is indeed the case: conceptually these phase 



(/) = exp [-t X 10 X (-(IE) + h + S,- 2IzS,)] . 



[35] 



(90° J,) (90° 5,) (-90° 2/,^,) 



[36] 
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shifts can be thought of as arising from the lack of a term in controlled phase 
shift gates, and thus have the desired properties. 



3.3 The controlled square- root of not 



Implementation of the controlled square-root of not gate is simple using the 
approach outlined in equations BBI and 133, with cj) = 7t/2. The controlled phase 
shift is simply implemented as 



(45°/,)(45°5,)(-45°2/,5, 



[37] 



Note that there is a close relationship between this gate and the spin-state se- 
lective excitation sequences (||^||^ which have been suggested as a method for 
simplifying E-COSY spectra. Clearly any other spin-state selective pulse can be 
created in much the same way. 



4 Three qubit gates 

A wide variety of three bit gates have been investigated, but we will confine our 
discussions to the Toffoli gate ([7^, or controUed-controUed-NOT. This takes the 
form 

/I 0\ 
1 
1 
1 
1 
1 
1 
VO 1 0/ 



- N - 



[38] 



and plays a central role in the theory of classical reversible computation as it 
can be shown to be universal (that is, any reversible classical logic circuit can be 
constructed entirely out of classical Toffoli gates). 

It might seem that this gate could be implemented using double-controlled 
phase shifts. 



- N - 



( 


» 










F 










h-1 




77 




h 





[39] 



but while the above circuit is indeed correct this approach is not practical. The 
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double-controlled phase shift matrix, tt can be decomposed as 



TT = exp 



Rz — Sz + 2IzRz + '2.IzSz + "^RzSz — AIzRzSz 



[40] 



but this decomposition cannot be used as a guide to implementing tt, as there 
is no NMR Hamiltonian directly corresponding to AIzRzSz- For the same reason 
it is not possible to directly implement a double-controlled square-root of not, 
that is a doubly spin state selective excitation sequence. 

It is, however, possible to implement these gates by using more complex 
networks of logic gates. Indeed it has been shown that the combination of a 
controUed-NOT gate and a set of general one bit gates is universal p^), so that 
any other gate can be constructed from them. This process is even simpler if 
the set of basic gates is slightly expanded; for example a Toffoli gate can be 
implemented using the network 



N 



N 



[41] 



N 



V 



-V- 



V 



It should be noted that while it will be difficult to construct a true Toffoli gate 
it is relatively simple to construct an approximate Toffoli gate, whose transforma- 



tion matrix has the same underlying form as equation |3^, but where the non-zero 
matrix elements are not all equal to unity. Such a gate was demonstrated early 
on (H), and can be used instead of a true gate in some situations where it is the 
last gate in the logic network, such as error correction 



5 Applying logic networks in NMR 



Having shown that NMR pulse sequences can be used to implement all the basic 
logic gates needed both for classical reversible computation and for quantum com- 
putation, it is instructive to consider whether any more complex logic networks 
might correspond to interesting NMR pulse sequences. One obvious candidate is 
the double-controlled square- root of not circuit, which corresponds to a doubly 
spin-state selective excitation sequence. More generally this family of gates in- 
volves the creation of an effective Hamiltonian containing a term proportional to 
AIzRzSz, which is not easily accessible by conventional means, and may be useful 
in the generation of multiple quantum coherences. 
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Another logic network which might prove useful is the swap network 




N 



- N 



N - 



'1 0^ 

10 

10 

,0 1, 



[42] 



which completely interchanges the states of the two spins involved. Clearly this is 
closely related to sequences like the double INEPT transfer step used in sensitiv- 



ity enhanced HSQC experiments ([^, \2li, . Unlike conventional heteronuclear 
transfer steps, however, this sequence preserves the states of both spins, by per- 
forming a complete swap. 

More speculatively it may be possible to use quantum error correction codes 
{2(i,3i to reduce the effects of spin-spin relaxation upon NMR spectra. Initial 
experiments in this direction ([^ suggest that error correction does work in 
NMR experiments, but that it is unlikely to have much practical significance; 
this assessment may, however, prove too pessimistic. 
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